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Abstract 

We apply density functional theory to study the influence of polydispersity 
on the stability of columnar, smectic and solid ordering in the solutions of rodlike 
macromolecules. For sufficiently large length polydispersity (standard deviation 
0" > 0.25) a direct first-order nematic-columnar transition is found, while for smaller 
0" there is a continuous nematic-smectic and first-order smectic-columnar transition. 
For increasing polydispersity the columnar structure is stabilized with respect to 
solid perturbations. The length distribution of macromolecules changes neither at 
the nematic-smectic nor at the nematic-columnar transition, but it does change at 
the smectic-columnar phase transition. We also study the phase behaviour of binary 
mixtures, in which the nematic-smectic transition is again found to be continuous. 
Demixing according to rod length in the smectic phase is always preempted by 
transitions to solid or columnar ordering. 

PACS numbers: 87.15. Da, 64.70.M, 61.30.Cz. 
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It has been known for some time that concentrated solutions of DNA-^, po- 
lypeptides^'^, polysaccharides^ and hairy rod polymers'^ form columnar phases. It 
also has been observed that DNA in bacteriophages and sperm nuclei of sepia, 
trout and salmon exhibit columnar ordering-^ . The activity of DNA (renaturation, 
transcription or replication) can be enhanced in the condensed phase^. Also, the 
condensed form of DNA can be used by nature to store genetic material in small 
volume and use it at the moment of cell cycle. Despite the accumulating body 
of experimental data, very little is known about the influence of various factors, 
such as attractive forces or polydispersity, on the stability of columnar ordering 
in macromolecular solutions. It is known that columnar ordering is preempted 
by smectic ordering in hard monodisperse rod systems^, but can be stabilized in 
binary mixtures of rods of different length^. However, since in actual solutions the 
systems are characterized by the continuous distribution of molecular length, true 
monodisperse or bidisperse systems are rare. Here we fill this apparent gap and 
study, within Density Functional Theory^'^ (DFT), the influence of polydispersity 
on columnar ordering. 

A polydisperse system with continuous distribution of molecular masses (or 
length as is the case here) can be regarded as a mixture of inflnitely many compo- 
nents. Thus, phase equilibria between two phases requires the equality of chemical 
potentials for molecules of all lengths, making the total number of equilibrium con- 
ditions inflnite. Previous studies of polydisperse systems have employed bifurcation 
analysis^° or expansion in the distribution function width-^-*^ (valid for sharp dis- 
tributions only), but, to date, no general approach is known. Here we establish 
the equilibrium conditions by a novel technique that does not involve expansion or 
assumption of the sharpness of the molecular length distribution. 

We pose the following questions: What is the minimal degree of polydispersity 
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necessary to stabilize the columnar phase? At what polydispersity is the smec- 
tic phase destabilized? Does polydispersity change the continuous nature of the 
nematic-smectic phase transition? Are there any smectic phases in which rods of 
different lengths are demixed (completely demixed, or ordered in layers with varying 
widths)? 

We consider a polydisperse system of hard, parallel, cylinders of diameter D 
interacting via the hard core repulsion potential. The free energy of the system as a 
functional of the number density po(r, L), for a given length LI (L is dimensionless) , 
is given by the following formula (the free energy is in units): 

F[p] = j dL j drpoir, L) (in (AVo(r, L)) - l) + j dL j drpo{r, L)^ {voPo{r, L)) . 

(1) 

where A is the De Broglie wavelength. The first term in Eq(l) is exact and represents 
the configurational entropy of the polydisperse system. The second term is the 
excess free energy determined by the interparticle interactions. Here is the excess 
free energy density of the homogeneous system and po{r, L) is the weighted density, 
defined by 

po(r, L) = j dL' j dv'w{v - r', L + L')po(r', L'), (2) 
where the weight function w{r — r' ^ L -\- L') is normalized according to 

J drw{r-r',L + L') = l (3) 

The weighted density (Eq(2)) represents the influence of the total density of particles 
on the density of rods of length L at point r. The normalized length distribution 
of the cylinders is given by: 

^(L) = l J drpoir, L), (4) 



where N is the number of particles in the system. 

Now we make the following approximations. First, for the weight function we 
choose the normalized Mayer function, i.e.: 

Q{{L + L')l/2-\z-z'\)e{D-\r^-r'^\) 
w{r-r,L + L)- ^^-^^ -^^ 

=hiL + L',\z-z'\)f2iD,\r^-r'^\), 
where is the Heaviside step function. This weight function implies that the 
free energy of a system with columnar ordering does not depend on the degree of 
polydispersity. Second, we employ a decoupling approximation, i.e., 

po{r,L) = p{r)g{L). (6) 

For the columnar phase Eq(6) is actually exact. For the smectic phase the ap- 
proximation is expected to break down only if the width of the layers varies. As 
we show below (for binary mixtures) however, this does not happen. The density 
distribution p(r) is approximated by a Gaussian function centered at the sites of 
the Bravais lattice characteristic for the given phase. It follows from the studies of 
hard spheres fluids®' -^^"-^^ that this is a very good approximation. We assume that 
g{L) is given by the Gaussian distribution, characterized by the mean value Lq, and 
the standard deviation a. It turns out that the mean length Lq scales out. Finally, 
the excess free energy density for the homogeneous system is approximated by the 
Carnaham-Starling equation as in Ref.6, i.e., 

*(r7) = ifclM (7) 
(l-r/)2' 

where rj — pvo, p is the average number density and vq — {1/4:)'kD^ILq is the mean 
volume of the cylinders. 

The phase diagram obtained from the above outline is shown in Fig.l. It 
encompasses the nematic, smectic, columnar and solid phases. For a > 0.25 there 



is a direct first order phase transition to the columnar phase. The chemical potential 
as a function of L for the nematic and columnar phases can be written in the form 



(8) 



It is central to our approach that for at least one of the phases in equilibrium g{L) 
may be calculated as a function of //(L), as is the case for Eq(8). We make use of 
this in Eq(13). In the nematic phase po is given by 



The equilibrium density distribution in the columnar phase p'^°^{r±) is a sum of 
Gaussian functions centered at the sites of a hexagonal lattice. The lattice constant 
and the width of the Gaussian peaks are obtained from minimization of the func- 
tional (Eq(l)) with respect to these variables. It can be seen directly, that the distri- 
bution function does not change at the nematic-columnar coexistence. The densities 
at coexistence, normalized by the density at close packing (rjcp = 7r/2-\/3 ~ 0.907), 
arc determined as follows: rf^'^^ /r]cp — 0.36 and r7*^°V^cp — 0.43, independent of the 
polydispersity of the system. 

For standard deviations a < 0.25 we find a continuous transition from the ne- 
matic to the smectic phase. The length distribution does not change the continuous 



(9) 



and in the columnar phase by 




(10) 
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nature of the transition. This result is in agreement with computer simulations 
by Stroobants^, who found a continuous nematic-smectic transition for a binary 
mixture of long and short spherocylinders. This result is not at all obvious. In 
principle one might expect that the nematic-smectic transition could be accompa- 
nied by the separation of rods, similar to the isotropic-nematic transition. In the 
latter case the longer rods are more abundant in the nematic than in the isotropic 
phase^^'^-*^'^^'^^. By analogy, we could expect the length distribution to narrow at 
the nematic-smectic transition and consequently to change the continuous nature of 
the transition. The transition has been studied as follows. First we have calculated 
the chemical potential of the smectic phase. It reads: 

//^-(L)=ln((7«-(L)) + 

/ dzp^'^iz) In (A V^(-s)) / d2p'"'(^)^'(i'op'™(-s, L)) 

J dzp'-'^{z) / dzp'-'^iz) (11) 

JdzJ dz' J dL'fi{L + L', \z - z'\)p'''^{z')vop^"'{z)g^"'{L')'^'{vop'"^{z' , L')) 

J dzp'''^{z) 

where f\ is defined in Eq(5). Then we have equated nematic and smectic chemical 
potentials: 

^nem^^) = ^«'"(L). (12) 

It follows immediately, that for the known distribution function g^^{L), the distri- 
bution function g^^^{L) is trivially determined at coexistence by Eq(8), i.e., 

^--(L) = exp (/x^-(L) - /.--) (13) 

The second equilibrium condition, i.e., the equality of pressure, combined with the 
normalization condition J dLg^^^{L) = 1 determines the coexisting densities. We 
have found further that for all degrees of polydispersity a the nematic-smectic phase 
transition is continuous. 
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The same procedure has been apphed to the columnar-smectic phase transi- 
tion. In this case we have assumed a given length distribution function in the 
smectic phase and from the equahty of chemical potentials (Eqs(8,10,ll)) we could 
determine the distribution function in the columnar phase at coexistence with the 
smectic. This transition is first order, so the average volume fraction rj jumps at 
the transition. The distribution function in the smectic is again assumed to be 
Gaussian. We characterize the distribution function in the columnar by its mean 
value and standard deviation in order to compare it to the distribution function in 
the smectic. It turns out that the change in the mean length is negligible. The 
standard deviation of the distributions is larger in the columnar than in the smec- 
tic. This is the expected result, since the lamellar ordering in the smectic favours a 
sharp distribution whereas the columnar structure does not. For example, as seen 
in Fig.l, a smectic phase with a polydispersity of cr = 0.15 at a packing-fraction 
of r]/rjcp = 0.45 is in equilibrium with a columnar structure of a polydispersity of 
(J = 0.17 at a packing-fraction of rj/rjcp = 0.50. 

The dashed lines shown in Fig.l represent the instability of the smectic and 
columnar phases with respect to perturbations to a hexagonal solid. The phase 
characterized by the density distribution p*(r) is unstable with respect to the per- 
turbation (5p*^*'-^)(r) of the symmetry of the / phase if the following condition holds: 



We have assumed that the perturbations can be expressed in the factorized form: 



where dp-^ (r) describes the onset of the ordering specific for the phase / and absent 
in the phase i. For analysis of the columnar solid bifurcation we take 




(14) 



(15) 



5pf{r) = cos (kz) 



(16) 
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while for the smectic sohd bifurcations we assume 

3 

5p^(r) = J]cos(k,r^), (17) 

n=l 

where ki = {l,l/V3)k, /ca = {-l,l/V3)k and ks = (0,2/\/3)/c are the vectors 
spanning the first shell in the reciprocal space for the regular hexagonal lattice. 

For large polydispersity the columnar phase stabilizes with respect to solid per- 
turbations. We expect this, since in the system of rods with continuous distribution 
of length, the particles do not fit well into the 3D structure involving ordering along 
the long axis of rods. This transition is expected to be first order, with the distri- 
bution function more strongly peaked in the solid. To quantify this in detail, it is 
necessary to examine the columnar-solid coexistence. The smectic phase becomes 
slightly destabilized with respect to solid perturbations if the polydispersity is in- 
creased. This can be understood as an indication that order in the z-direction is 
least favoured, so that even 3-dimensional ordering is preferred. 

To complete this study, we investigate a binary mixture of rods of two different 
lengths. The aim is to find out, under what conditions we have to expect any form of 
demixing. Since a bidisperse system is more likely to segregate into its components 
than a polydisperse one, we consider this to be a stronger criterion than looking for 
an instability of a polydisperse system. We have already found, that no demixing 
occurs at the nematic-smectic transition in the polydisperse system. This result 
remains also valid in the binary mixture. 

The functional given by Eq(l) for the general case of a polydisperse system 
reduces to the binary mixture case for 

P(r, L) = pi(r)5L,L, + p2{r)dL,L2, (18) 

where 5ij is the Kronecker delta function. Using this density distribution, we are 
able to calculate the free energy of the nematic, smectic and columnar phase as 
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well as the instabilities to solid ordering as before. Fig. 2 shows the results for 
different ratios of the lengths of the rods. The phase-diagram is calculated at the 
equivalence point, where the partial volume fractions of the two components are the 
same. Although the results are not completely comparable to those of Stroobants^, 
who studied spherocylinders rather than cylinders, qualitative agreement can still 
be seen easily. As in the study of Stroobants, we observe the destabilization of the 
smectic order compared to the nematic, and stabilization of the columnar order 
with an increasing length-ratio of the two components. The qualitative similarity 
of the phase diagrams for the polydisperse and bidisperse systems is remarkable 
(compare Fig.l and Fig.2). This is a quite unexpected result and indeed may have 
interesting implications for the future study of such systems. 

We now consider the possibility of demixing within the smectically ordered 
system. One might, for example, expect that alternating layers of different widths 
would be formed in the bidisperse system. We look only for a separation of the 
complete system into one phase consisting mainly of long rods and another phase 
consisting mainly of short rods, since it can be shown, that the free energy of such a 
system is comparable to that of smectics with alternating layers. In order to study 
the stability we check the partial derivatives 

= (19) 

where pi are the partial densities of the two components. It turns out that the ma- 
trix (Mij) remains positive definite for all length-ratios, compositions and packing- 
fractions, for which the smectic phase is stable. We also check for the possible coex- 
istence between two smectic phases with widely varying composition, by equating 
the chemical potentials of both components as well as the pressure. Again we find 
no such coexistence in the smectic phase. This leads to the conclusion that demixed 

9 



smectics or smectics with varying lamella-widths are not found in binary systems 
under the present conditions. It follows that demixed smectics are unlikely to occur 
in multi-component systems. Demixing may well occur in the solid, however. 

A physical realization of our model is a dense polydisperse suspension of elon- 
gated colloidal particles stabilized against irreversible aggregation by surface-grafted 
polymer layers. In a good solvent the interactions between such colloidal particles 
can be approximated by hard core interactions^°'^^. 

In general, however, macromolecules in solutions interact via van der Waals 
forces, and if they carry surface charges, also via screened Coulombic repulsion^°'^^. 
The influence of these complex interactions on the phase diagram is not known in 
general. 

Different problems arise when the macromolecules are not rigid but semi- 
flexible^^'^^. It should be noted as well that in self assembling systems, when the 
length distribution is given by the thermodynamic conditions in a given phase^^'^^, 
the equality of monomer chemical potential is sufficient to set the coexistence con- 
ditions between two phases. 

Summarizing; by means of density-functional theory, we have studied in de- 
tail the influence of polydispersity on phase equilibria in oriented rodlike macro- 
molecular systems. For a Gaussian distribution of rods length we find three poly- 
dispersity regimes. For 0.15 > cr > we have the following sequence of phase 
transitions nematic-smectic-solid; for 0.25 > u > 0.15 we find the nematic-smectic- 
columnar and solid phases; for 0.27 > a > 0.25 we find direct nematic-columnar and 
columnar-solid transitions and for a > 0.27 there are only nematic and columnar 
phases. We also find that polydispersity does not affect the continuous nature of the 
nematic-smectic phase transition. We do not find any evidence for demixed smectic 
phases. Our novel method for the study of polydispersity can be easily applied also 
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to the isotropic- nematic phase transition^^. 
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Figure Captions 

Fig.l Phase diagram for a polydisperse system of parallel rods interacting via hard 
core repulsive forces. The polydispersity of the rod lengths is modelled by a 
Gaussian distribution of standard deviation a. Squares denote coexistence of 
phases, triangles a second order phase transition, and crosses represent insta- 
bilities. 

Fig.2 Phase diagram for a two component system of length ratio L1/L2. Note the 
remarkable qualitative similarity to Fig.l. 
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